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A model for correlated electrons in a lattice with local additional spin-1 degrees 
of freedom inducing constrained hopping, is studied both in the low density limit 
and at quarter filling. We show that in both ID and 2D two particles form a bound 
state even in presence of a repulsive U < U c - A picture of a dilute Bose gas, leading 
to off-diagonal long range order (LRO) in 2D (quasi-LRO in ID), is supported by 
quantitative calculations in ID which allow for a determination of the phase diagram. 
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Recent work on models with correlated electronic hopping, involving density-dependent hopping 
matrix elements, has pointed out that superconductivity can arise by purely off-diagonal processes 
in spite of repulsive short-range interactions. 

Correlated electronic hopping arises also when additional degrees of freedom participate in the 
"local physics" at each site, inducing, for instance, occupation- dependent constraints. One clear 
example of the latter situation is realized in molecular crystals, where the additional degrees of 
freedom arise because of the coupling to intra-molecular phonons. pL[j| The basic physics is that, 
under appropriate conditions, molecules with an odd number of electrons have an associated Berry 
phase and an orbital degeneracy, due to a dynamical Jahn- Teller effect, those with even number 
have none. Q The general question is whether such mechanisms will also lead to superconductivity 
even in the presence of repulsive diagonal interactions. 

The simplest model capturing some of the physics of electronic hopping in presence of additional 
local degrees of freedom is a Hubbard-type model with an extra spin-1 at each site. || The constraint 
is that the allowed spin-1 states depend on the electron occupancy: [|| an empty or doubly occupied 
site must have S z = 0, whereas a site occupied by a single electron (either up or down) has an 
additional twofold (orbital) degeneracy, represented by S z = ±1, 

n r = 0, 2 -> S z = , n r = 1 S z = ±1 . (1) 

The hamiltonian H is written as follows: 

H = lEE (4<r<V* + H - c -) ( S r S r> + H.c.) + C/^n rT n ri , (2) 

<rr'> a r 

where the Sy s are spin-1 ladder operators at each site, and the remaining notation is completely 
standard. It is worth stressing that while H, and in particular the hopping term, conserves the 
constraint in Eq. [l], the model is still highly non-trivial, even for [7 = 0. 

In this paper we will present some results concerning the few particle problem, the low density 
region, and the quarter-filling case, which shed light on the basic features of the phase diagram of 
the model. We will show that in both ID and 2D the additional spin-1 degrees of freedom lead to 
a two-electron bound state even in presence of a repulsive U < U c , Moreover, the model does not 
show tendency to phase separation and a picture of a dilute Bose gas, leading to off-diagonal quasi- 
long-range-order (LRO) in ID (genuine LRO in 2D) is supported at small density by analytical 
calculations. At quarter-filling, a clear superconducting regime is found to survive up to a positive 
U~t. 

Consider the two-particle problem first. A state in the two-particle Hilbert space with total z- 
component of the spin M^ OT = (for both the electron spin and the spin-1 states) can be written 
as |*) = Ii„,^ + -(r,r')S+5; + ip-+(r, r')S~ S+jc^c^jjO), where the vacuum |0) is the state 
without fcrmions and with S z = at each site. In writing we have taken into account the two 

possibilities of associating a S z = ±1 spin state to the up and down electrons: i/'h is the amplitude 

for having S z = +1 associated to the "["-electron (and S z — — 1 to the |-electron), while ip |_ is the 

amplitude for other possible choice. The Schrodinger equation for ip_\ (r, r') is easily shown to be 

Etp + - (r, r') = -t^2[ip + -{r + a,r')+ip+- (r, r' + a)} + U 5 r y tp+- (r, r') 

a 

-t(J2Sr+a,r>) [ r/>-+(r, r) + V-+(r', r') ] , (3) 

a 

where a denotes a nearest neighbor vector (a = ±1 in ID). A similar equation is obtained for 
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ip r') by just exchanging ijj |_ and ip^ everywhere. The last term in Eq. || is crucial to 

the whole story, and deserves a few comments. When the two electrons are far enough in the 
otherwise empty lattice, the hamiltonian H simply allows the hopping to a nearest neighbor site 
of the "composite" object formed by an electron and the associated spin-1 state (first term in Eq. 

. Things are more subtle when two electrons come to the same site r. In such a case, from a 
doubly occupied site with S* = one can reach, upon hopping, two possible final states: either 
each electron keeps its own spin-1 state or the spin-1 states associated to the two electrons are 
exchanged. It is precisely this second possibility of exchanging spin-1 states that is responsible for 
the presence of tp |_ in the equation for ijj^ and vice- versa (last term in Eq. ^|). 

The Schrodinger equation is easily solved in momentum space, where it reduces to a 2 x 2 matrix 
problem. For total momentum P = 0, the set of solutions among which the ground state is found 
satisfies the equation 

— y^— = — (4) 

k 



where E is the energy eigenvalue and e& is the tight-binding dispersion of the free-electron problem 
(efe = — 2t J2 a cosk a ). In the ordinary Hubbard case, the right hand side of Eq. |] would simply read 
l/U. H A graphical analysis of Eq. ^ readily shows that a bound state solution is present even for 
U > up to U c — ADt in D < 2. In D > 3 a finite attractive U is needed to produce a bound state. 

The bound state solution can be worked out analytically in ID. For general values of the total 
momentum P and for U < U c = At cos(P/2), the ground state energy is given by 

E P = -At cos (P/2) (Z 2 P + 1)/(2Z P ) , (5) 



with Z P = {-(U/At) + [(U/At) 2 + 3 cos 2 (P/2)} 1 / 2 }/ cos (P/2), and the corresponding ground state 

wave-function is ijj^ = ip |_ = e lP ^ r+r )/ 2 [e~ K l r ~ r I — (\/2)8 r _ r i] with k = In Zp. At larger values of 

U, no bound state is present and the energy spectrum is continuous in the infinite lattice. There is 
still, however, an antibound state. The two-particle solution cannot be generalized to an arbitrary 
number of particles by Bethe Ansatz because the corresponding scattering matrix does not satisfy 
the Yang-Baxter relations. 

The form of the ground-state wave function naturally provides a picture of bound pairs approx- 
imately localized on adjacent lattice sites, thereby forming dimers. Remarkably, the rather strong 
attraction responsible for this binding is generated by the kinetic term alone via the presence of 
the additional degrees of freedom. The critical value of the Hubbard repulsion U c — At is consid- 
erably larger than the ground state binding energy at U = 0: Eb — At(Zo — 1) 2 /(2Zq) 0.618 1 
showing that this kind of pairing mechanism is rather insensitive to the presence of on site Coulomb 
repulsion. This is explained by its "kinetic" origin which delocalizes the pair on neighboring sites. 
The same feature is also present in 2D where U c = 8t: The enhancement is due to the the larger 
coordination of the 2D lattice which provides an even more efficient derealization of the electron 
pair. 

This interpretation of the two-particle ground state in terms of dimers leads to a simple picture, 
both in ID and 2D, of the low density region of the model (||), a picture which we will explicitly 
verify in ID: For U < U c , at low densities, the system behaves as a weakly interacting, dilute gas 
of dimers, which follow boson statistics and may be thought of as bosons with an extended core. 
In ID, we can estimate the dimer effective mass from our two-particle calculation, Eq. pi which 
at small momentum P gives Ep ~ Eq + t c ^P 2 with t e g = t [3 + (U/At) 2 ~\ . The effective mass 
has a smooth dependence on U varying between y/3 at U — and 2 at U = At. Dilute dimers are 
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therefore rather mobile and a superfluid ground state must be expected at zero temperature. In 
ID, of course, off-diagonal LRO cannot occur and only a long range power-law decay of the dimer 
density matrix is possible, while in 2D a genuine Bose condensate will form. In terms of the original 
electrons this implies a standard strong coupling BCS superconducting ground state with localized 
Cooper pairs. A similar scenario has also been proposed in the framework of the one dimensional 
t — J model where bound pairs are formed at low density 0] and 2 < J/t < 2.95. In that case, 
however, the model is unstable to phase separation, which in fact occurs massively at larger values 
of J/t. 

To check whether the superfluid dimer picture is correct and that no phase separation occurs 
at low density, we have carried out a numerical investigation of the four particle problem by exact 
diagonalization of the model (0) in lattices with up to 24 sites. The results for the ground state 
energy and for the density-density correlation function at U = are reported in Fig. 1. The four 
electrons are sharply localized in two pairs placed at the maximum separation. The origin of this 
repulsion is purely kinetic and is also reproduced by the ground state wave function of two hard 
core bosons in ID, tp(R) °c | sin (7ri?/L)|. The scaling of the energy of four electrons clearly shows 
that the limiting value is just twice the pair energy E = — 8£/v3 with l/L 2 corrections. Free "hard 
core bosons" with the appropriate effective dispersion t c g = t [3 + (t//4t) 2 ] , give the asymptotic 
result L 2 (E — 2Eq) — > 27r 2 i e ff marked with a full circle in Fig. 1, and clearly compatible with the 
data. 

A further step in the quantitative characterization of the low density region of the model can 
be made by using standard field theoretical methods for ID systems. Q The analysis of the two 
and four-particle problem shows that for U < U c a spin gap is present in the excitation spectrum 
of the model while the charge (superfluid) degrees of freedom remain gapless, giving origin to the 
power-law scaling of the ground state energy shown in Fig. 1. Above U c a two-electron bound state 
does not form and the triplet spectrum looks gapless. This is actually not the case at finite density, 
as can be inferred from the analysis of the U — > oo limit of the model. |Q The strong coupling 
expansion of the model in Eq. ^| can be obtained closely following the analogous calculation in the 
Hubbard model. Q The resulting hamiltonian, to 0(t 2 /U), is defined in the subspace of singly 
occupied sites where electrons are characterized by two internal spin-1/2 degrees of freedom: The 
usual spin (a) and a pseudospin (r). The latter takes into account the two possible states of the 
original spin-1 allowed by the constraint (0). The resulting hamiltonian contains a free hopping 
term (with coupling — t) and a spin (and pseudospin) dependent part (with coupling oc t 2 /U). In 
ID, analogously to the Hubbard model case, (l(J the wave function factorizes, at U — > oo, for all 
densities: The position of the electrons is determined by a spinless fermion wave function while the 
spin (and pseudospin) ordering is governed by the effective spin hamiltonian 

H ST = -J cS ■ B i - R ' ^ - V 4 ] » (6) 

<i,j> 

with J e g depending on the electron density p = N/L as J e g = (8t 2 /U)p[l— sin(27rp)/27rp]. Numerical 
as well as variational calculations have shown that the hamiltonian in Eq. || is characterized by a 
spin gap of order J c g. R] Fig. 2 (inset) shows the finite size spin gap at filling p = 1/2 for several 
values of U . A finite spin gap is quite clear from the L — > oo extrapolations. We also find that the 
gap persists to large U and correctly scales with 1/U. Similar results are also obtained at half-filling, 
where Eq. || is the large U effective hamiltonian. We conclude that our model (Eq. |^), at arbitrary 
(finite) density, has a spin gap all the way to U — > oo. In the zero-density limit p — > 0, however, 
the spin gap vanishes very fast, like p 3 , whence the gapless appearance of our previous two-particle 
results. 

The spin (and pseudospin) degrees of freedom are gapped at all U and only the charge degrees 
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of freedom remain gapless. The universality class of the hamiltonian in Eq. is then given by the 
Luther-Emery model, [0 and the long range decay of the correlation functions is characterized by 
power-laws which can be obtained by the knowledge of the single exponent K p . |(| This exponent 
is also related by the Haldane-Schulz equation to the compressibility of the model: 

d 2 E _ 7T u p 



The factor 4 at right hand side takes into account the presence of an additional (pseudospin) degree 
of freedom which can be attributed to the electron in the present model, due to the additional 
degeneracy S z — ±1 in Eq. [I]. Correspondingly, N includes both spin and pseudospin degrees of 
freedom. As usual, u p represents the charge velocity, i.e., the velocity associated with the branch 
of gapless excitations. In order to use Eq. for obtaining the low density limit of K p , we should 
calculate both the ground state energy and the charge velocity at small but finite density. However, 
assuming a smooth zero density limit, it is possible to get the same information from the two 
and four-particle solution. In fact, at fixed JV, and L — > oo, the predicted scaling of the energy 
is E(N, L) = Ne + aN 3 /L 2 - (3N/L 2 + 0(L^ 3 ). Analogously, the large L (i.e., low density) 
limit of the charge velocity should scale as u p — jN/L, where a, (3 and 7 are model-dependent 
constants. By means of these asymptotic expressions together with Eq. (?], we get the formal result 
for K p — 7r7/24a. The central charge c jnj of the model can be also expressed in terms of these 
constants, c = 6[3/irj. This will provide a consistency check. At U > U c — At, the two-particle 
analytical result is sufficient to extract the values of the required constants, which turn out to be 
independent of U: a = 7r 2 /24, f3 = 7r 2 /6 and 7 = ir. This gives c — 1, as expected from the 
universality class of the model, and K p = 1 . In order to extract the correct large L scaling of the 
energy in the regime U < U c , the four-particle numerical result is required. By using the estimate 
of Fig. 1, the parameters governing the asymptotic low density behavior become a = t ff7r 2 /192, 
= £cff7r 2 /12 and 7 = f e ff7r/2, which give c = 1 and K p — 4 independent of U < At. 

The exponents which characterize the decay of the electronic correlation functions are obtained 
in terms of the parameter K p by bosonization under the assumption that only the charge sector 
remains gapless. || In particular, the density response function at 2Uf decays as x ~ Kp / 2 while 
the Akp component behaves as x~ 2K, >. In the derivation of these power-laws, we must recall that 
eight bosonic fields contribute to the physical density operators, and this reduces the exponents 
by a factor 2 with respect to the usual Hubbard model case where only four fields occur (right 
and left movers with two possible values of the spin). At the Luther-Emery fixed point also the 
superconductive correlations have power-law behavior. By examining the corresponding operator, 
we get the asymptotic behavior of the s-wave superconductive correlation as A s (a;) ~ x~ 2 / K >>. At 
low densities and for U < At we have K p = 4: the most relevant correlations are the superconductive 
ones, decaying as 1/^/x. Remarkably, the 1/y/x behavior is identical with that of a dilute hard core 
boson gas, as expected from the physical picture previously discussed. Instead, for U > At and 
p — * 0, we have K p = 1: the correlation with the slowest decay is the density response function 
at 2kp , with an inverse square root behavior, implying a divergence in the CDW susceptibility. 
Strictly at U = 00 the 2fc^ response function has vanishing amplitude (like in the Hubbard model) 
and therefore the corresponding singularity is absent in the exact solution previously discussed: we 
expect that, at finite U , this additional feature may be detected. 

We have also determined K p numerically from the charge velocity u p and from the Drude peak 
D = UpKp/(2ir). H Fig. 2 shows the results obtained at quarter-filling (p = 1/2). For U/t < 1 
the L — + 00 extrapolation clearly gives K p > 2, i.e., predominantly superconducting correlations. 
On the contrary, for U/t > 1 we have K p < 2, i.e., predominant CDW correlations. The crossover 
between the two regimes is found to take place close to U w t at this density. Exactly at half-filling 
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(p = 1) a transition is found at U c ~ between a superconductor ([/ < U c ) and a Mott insulator 
fully gapped in both spin and charge sector. j5j 

The physics of the model in one dimension is now rather clear. A transition occurs at U c (p) (with 
U c (0) = 4f, U c (l/2) w i, and I7 C (1) ~ 0) between a superconductor (U < U c ) and a sliding CDW 
conductor (p ^ 1) or Mott insulator (p = 1). The physics at small U, leading to superconducting 
correlations, is rather robust and can be generalized to 2D where true LRO will occur at zero 
temperature. On the other hand, the CDW state, perhaps akin to a spin Peierls state, should get 
weaker and disappear in higher dimensions. The large U spin gap is also, probably, a II? peculiarity. 

By extrapolating these results, we can surmise that a correlated hopping model such as Eq. 2 
can be expected to give rise to superconductivity in 3D. A pairing mechanisms based on correlated 
hopping is not easily destroyed by a repulsive U, it is more effective at low carrier density, and is 
apparently immune from the polaron self-trapping, which depresses T c in strongly coupled electron- 
phonon systems. More work is now being devoted to investigate the relevance of this type of model 
to molecular superconductors such as organic ID metals, Chevrel phases, and metal fullcridcs. 
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Figure Captions 

Figure 1 Size scaling of the ground state energy of the hamiltonian in Eq. 2 with U — and four 
electrons. AE = £?l(4) — 2E O0 (2) where El{N) is the ground state energy of N electrons 
in a L-site ring in units of t. The full circle is the exact result for two hard-core bosons in 
the infinite lattice with t e s = t/V%, fully compatible with the L — > oo extrapolation of the 
data. Inset: the density-density correlation function for four electrons in a L = 20-site ring at 
U = 0. 

Figure 2 Size scaling for the exponent K p , obtained from finite size values of the Drude peak D 
and charge velocity u p at filling p = 1/2, for U/t = 0,1,2, and 4. Inset: the spin gap for 
the same parameters. The continuous lines represent fits of the form AE = a + b/L + c/L 2 
and K p = d + e/L 2 . Note the transition from superconductive (K p > 2) to CDW correlations 
(K p < 2) for U/t w 1. The spin gap is always finite. 
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